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Abstract 

We propose a notion of operator monotonicity for functions of several vari- 
ables, which extends the well known notion of operator monotonicity for func- 
tions of only one variable. The notion is chosen such that a fundamental rela- 
tionship between operator convexity and operator monotonicity for functions of 
one variable is extended also to functions of several variables. 



1 Introduction and main result 

The notion of operator convexity for functions of several variables has been extensively 
studied in the literature. The first step is to define the functional calculus for functions 
of several variables. This can be done in the following way: 

Let 1% , . . . , Ik be real intervals and let / : 1\ x ■ ■ • x ly. — > R be a Borel measurable 
and essentially bounded function. Let x = (xi, . . . , Xf.) be a fc-tuple of bounded 
self-adjoint operators on Hilbert spaces Hi, . . . , Hk such that the spectrum of x% is 
contained in Ii for i = 1, . . . , k. We say that such a fc-tuple is in the domain of /. If 



Xi = J Xi Ei(dXi) i = 1, 
is the spectral decomposition of Xi, we define 



(1) /0)= / f{Xi,...,\ k )E 1 {dXi)»-"»E k (dX k ) 

Jl 1 x-xl k 

as a bounded self-adjoint operator on Hi ® •• • ® Hk, cf. [|[ [j], g]. If the Hilbert 
spaces are of finite dimension, then the above integrals become finite sums, and we 
may consider the functional calculus for arbitrary real functions. This construction 
extends the definition of Koranyi Q for functions of two variables and have the 
property that 

f(xi, ...,X k ) = fl{xi) <g> • • • <g> /fc(Xfc), 

whenever / can be separated as a product f{t\, ...,%) = /i(ii) • • • fk(tk) of fc func- 
tions each depending on only one variable. 

Remark 1.1 One might consider the functional calculus only for commuting opera- 
tors Xi, . . . , Xk on a single Hilbert space H and define 

fcom(xi, Xk) = I f(X 1 ,...,X k )dE(Xi,...,X k ) 
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as an operator on H, where E is the product measure of the commuting spectral 
measures associated with each of the operators. This approach was suggested by 
Pedersen and Lieb in [l^]. Our definition in equation (|l|) can then be written as 

f(xi, ...,x k ) = f com (x 1 <g> 1 <g) • • • <£> 1, . . . , 1 ® • • • ® 1 ® x k ) 

for arbitrary non-commuting operators x±, . . . ,x k on ii. If however the operators 
xi, . . . ,x k do commute, then there is a self-adjoint projection P on H (g) • • • (g) if 
with range isomorphic to ii such that f C om{x\, . . . , = Pf(x±, . . . , x k )P. The two 
approaches are thus essentially equivalent. 

Once the functional calculus is defined, we say that a function / : I\ x • • • x — > R 
is operator convex, if / is continuous and the operator inequality 

f(Xx + (1 - A)y) < A/(a;) + (1 - X)f(y) VA e [0, 1] 

holds for all fc-tuples of self-adjoint operators x — (x\, . . . , x k ) and y = (yi, . . . , y k ) in 
the domain of / acting on any Hilbert spaces Hi, ... , H k . The definition is meaningful 
since also the fc-tuple Air + (1 — X)y is in the domain of /. We say that / is matrix 
convex of order (rti, . . . , n k ), if the operator inequality holds for operators on Hilbert 
spaces of finite dimensions (m, . . . , n k ). 

The aim of this paper is to define the notion of an operator monotone function also 
for functions of several variables. The definition should, when restricted to functions 
of only one variable, be a simple reformulation of the ordinary condition for such 
functions. We also want the following theorem to be valid. 

Theorem 1.2 Let f : [0, a\[x • ■ • x [0, a k [— > R be a continuous real function. The 
following statements are equivalent: 

(i) f is operator convex, and f(ri, . . . , r k ) < if ri — for some i — 1, . . . , k. 

(ii) The function g :]0, a\[x ■ ■ ■ x]Q, a k [—* H defined by setting 

g(ri, ...,r k )= rV 1 ■ ■■r' k 1 f(r 1 , . . . ,r k ) 

is operator monotone. 

The theorem above is known to be valid for functions of one variable 0, 2.4 
Theorem] , and the extension to functions of several variables seems to be very natural. 
Our notion of operator monotonicity for functions of several variables is ultimately 



given in Definition 2.14, but it depends on intermediary notions and results given in 



Definition 2.1, Definition 2.2, Definition and Corollary 2.13. 

Before proceeding with this programme, we shall briefly discuss other possible 
definitions of operator monotonicity for functions of several variables, which we ulti- 
mately have rejected. 

Proposition 1.3 Let f be a non-negative continuous function of k variables defined 
in the first quadrant [0, oo[x • • • x [0, oo[. If f is matrix concave of order (n\, . . . , n k ), 
then 

0<Xi< yi i=l,...,k => f(x)<f{y) 

for arbitrary k-tuples of positive semi-definite matrices x = (x\, . . . , x k ) and y = 
{yi,...,y k ) of order (m, . . . ,n k ). 
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Proof: Let the appropriate fc-tuples of matrices be chosen and take A E [0, 1[. We 
set Zi = A(l — \)~ 1 (y i — x^ and notice that 

Xyi = Xxi + (1 — \)zi and Zi > 

for i = 1, . . . , k. Since / is matrix concave and non- negative we obtain 

/(Ay) > A/(:r) + (1 - X)f(z) > Xf(x) 

where z — (z 1; . . . , z^). The result now follows by letting A tend to one. QED 

The converse is not true. The function of two variables f(r±, r 2 ) — r\r 2 is indeed 
matrix increasing of any order in the sense that 

f(x 1 ,x 2 ) = xi <8 x 2 < yi <E> y 2 = / (3/1,2/2) 

for < x\ < yi and < x 2 < y 2 , but it is not even concave as a real function. How- 
ever, the situation is quite different for functions of only one variable. Mathias [Tl|| 
showed that a function, defined on the positive real half-line and matrix monotone of 
order n, is matrix concave of order [n/2]. It follows from ||, [7), although not stated 
explicitely, that a function, defined on the real positive half-line and matrix monotone 
of order An, is matrix concave of order n. If we relax Mathias' result only very slightly 
and are satisfied with proving that a function / : [0, 00 [— > R, matrix monotone of 
order 2n, is matrix concave of order n, then the following very simple argument will 
do. Let x\, x 2 be positive definite matrices of order n and notice || that to a given 
e > the inequality 

y* ( x i \ y = 1 ( x i + x 2 x 2 - £1 \ < / 2,- 1 (x 1 +x 2 )+e 
\ x 2 J 2 V x 2 - x 1 xi + x 2 J ~ \ A 

is valid for a sufficiently large A > 0, where 

\/2 / 1 -1 



V - 2 V 1 1 



is a unitary block matrix of order 2n x 2n. We then obtain 

1 / f(xi) + f(x 2 ) f(x 2 ) - f(x x ) \_y*( f(xi) 
2\ f(x 2 )-f( Xl ) f{x x ) + f(x 2 ) ) \ f(x 2 

^ \v\ < ( /(2- 1 (^i+^)+e) 

X2 y j - 1 /(A) 



f( X l) + f( X 2) < , ( Xl + x 2 



and consequently 



from which the statement follows by letting e tend to zero, since matrix monotone 
functions of order greater or equal to two are continuous (even continuously differen- 
tiable). 
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We shall finally mention that Koranyi and others have considered a notion of 
operator monotonicity for functions of two variables defined on I 2 where J =] — 1, 1[. 
The notion is closely connected to the theory of analytic functions of several variables, 
and in particular to a generalization of the Riesz-Herglotz formula |], [l^] . According 
to this theory the function 

g(ri ' r2)= (l + r?)(l + r 2 ) 



would be called operator monotone, but this is not consistent with Theorem L2 as 
the continuous function 

/(ri ' r2) = (l + nKl + r 2 ) ^ 2e[ °' 1[ 

is not operator convex. Koranyi's notion of operator monotononicity leads to no 
significant distinction between functions of one and two variables as does the theory 
of operator convex functions. 



2 Decompositions and monotonicity 

Definition 2.1 Let a; be a positive invertible operator acting on a Hilbert space 
H. We say that an /-tuple (y\,...,yi) of positive invertible operators on H is a 
decomposition of x (of length /) if 

(2) 2/i H \-yi = x. 

The /-tuple a — (ai, . . . , ai) defined by setting = x~ 1 / 2 y]^ 2 for i = 1, . . . , / is called 
the associated unitary row. 

We recall [0 that an /-tuple a = (ai, . . . , a;) of operators on a Hilbert space H is said 
to be a unitary row, if there exists a unitary operator U on the direct sum of / copies 
of H such that (ax, . . . , ai) is the first row in the / x / block matrix representation of 
U. The equation 

(3) aia\ + ■ • • + a;a ; * = 1 (the identity on H) 

is a necessary, but in general not sufficient condition for a = (a\, . . . ,ai) to be a 
unitary row. 

The row a — (ai, . . . , ai) associated with the decomposition of x in the definition 
above satisfy condition (||) since 

* i i * -1/2 -1/2 i i -1/2 -1/2 -, 

a\a l + ■ ■ ■ + aidi =x ' y\x ' -\ h x ' yix ' =1. 

Araki and the author proved that an /-tuple a — (ai, . . . , a;) satisfying condition (^) is 
a unitary row, if dimkcra^ = dimkera* for at least one i = 1, . . . , /. The condition is 
trivially satisfied in this case since the operators a%, . . . , ai are invertible. The /-tuple 
a = (ax, . . . ,ai) in Definition 2.1 is therefore indeed a unitary row. We notice that 
yi = a*xai for i = 1, . . . , I. 
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Definition 2.2 An index is a pair (l,j) of integers, where I > 2 and < j < I — 1. 

Definition 2.3 Let / :]0, ai[x ■ ■ ■ x]0, «fe[— * R be a real function. The constants 
ai, . . . , ate may be plus infinity. 

(i) We say that / is operator monotone of index if / is continuous and 

(*) dia,g(f(y tll ,...,y tkk )) < f(xi, . . . , Xk)L lk -i 

\ / \ t\—j (modi) 

for every fc-tuple x — [xi, . . . ,Xk) in the domain of / acting on any Hilbert 
spaces Hi , . . . , Hk and all decompositions 

Vu-\ Vyu=Xi i = l,...,k 

where L t k-i is the l k ~ 1 x l k ~ 1 block matrix with the unit operator on the tensor 
product Hi (g) • • • ® Hk in each entry. The index t is a multi-index of the form 
t = . . . , tk), where ti = 1, . . . , I for i = 1, . . . , k and weight \t\ = t\ + • • • 

[it) We say that / is matrix monotone of index (l,j) and order (ni, . . . , n/-), if the 
same inequalities (*) are satisfied for operators acting only on Hilbert spaces 
Hi, . . . , Hk of finite dimensions (ni, . . . , rik)- 

It is not difficult to establish that a continuous function is operator monotone 
of index [l,j), if and only if it is matrix monotone of index and all orders 

[ni, . . . , rife). The proof follows a suggestion by Lowner as reported by Bendat and 
Sherman Lemma 2.2] and can easily be adapted to the present situation. Fur- 
thermore, consider fc-tuples [mi, . . . , rrik) and (ni, . . . ,nk) such that m, < rij for 
i = 1, . . . , k. If a function is matrix convex of order (ni, . . . , n&) then it is also matrix 
convex of order (mi, . . . , mjt). Likewise, if a function is matrix monotone of index 
(Z, j) and order (m, . . . ,nk), then it is also matrix monotone of index [I, j) and order 
(mi, . . . ,m k ). 

Proposition 2.4 A continuous real junction f :]Q,a[— > R is operator monotone of 
any given index [l,j), if and only if it is operator monotone. Likewise is f matrix 
monotone of any given index [l,j) and order n, if and only if it is matrix monotone 
of order n. 

Proof: If we set k = 1, the inequality (*) reads 

/(Mi)</(*i) for j =0 

f{y 3 i)<f{xi) for j =l,...,l-l 

where yn + ■ ■ ■ + yn — xi is a decomposition of xi. These inequalities are trivially 
satisfied if / is operator monotone. If on the other hand one of the above inequalities 
are satisfied for a given index [l,j) and all decompositions of any xi in the domain 
of /, then / is operator monotone. The same reasoning applies to matrix monotone 
functions. QED 
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Proposition 2.5 Let f :]0,ai[x ■ • • x]0, ctk[— ► R be a continuous function and con- 
sider for i = 1, . . . , k the function of one variable 



9i(n) = f(r%, ■ 



,r k ) 



obtained from f by keeping all variables fixed except the ith variable. If f is operator 
monotone of some index then gi is operator monotone. Likewise, if f is matrix 

monotone of some index and order (m, . . . , n^), then gi is matrix monotone of 
order ni . 

Proof: Let / be operator monotone (or matrix monotone) of some index and 
assume i = 1. We choose operators y < x in the domain of g\. For some sufficiently 
small e > we set 



2/ml 



y 



e+(x-y)/(l-l) 
and for p = 3, . . . , k 



m = j + 1 



Urap 



and y m2 



m = I 
m =/= I. 



r-2 

e 



m = I — 1 
m / I - 1, 



We thus have the decompositions yn + ■ ■ ■ + yn = x + (I — l)e and 
yi P -\ 1- yi P = r p + {I - l)e p = 2,...,l. 

By only considering the index t = (j + 1, 1 — 1, 1, . . . , I) with length \t\ = j (modi) in 
(*), we obtain the inequality 

f(y, r2,-..,r p )<f(x + (l- l)e, r 2 + (I - l)e, . . . , n + {I - l)e) 

from which the inequality gi(y) < gi(x) is derived by letting e tend to zero. QED 



To further investigate the content of Definition 2.3 we set k = 2 and / = 2. The 
inequality (*) then reads 



and 



f(yn,yu) o 

f(V21,V22) 



f(yn,y22) o 

/(j/2i,J/ia) 



< 



< 



f(x 1 ,x 2 ) f{xi,x 2 ) 

f(xi,x 2 ) f(xi,x 2 ) 

f(xi,x 2 ) f(x 1 ,x 2 ) 

f(xi,X 2 ) f(x!,X 2 ) 







3 = 1 



for decompositions yn +y 2 i = x% and yi 2 +y 22 = x 2 . This is so because the solutions 
to the equation \t\ = t\ + t 2 = j (mod2) are the multi-indices (1, 1), (2, 2) for j = 
and (1, 2), (2, 1) for j = 1. If we set fc = 2 and 2 = 3 the inequality (*) reads 



/(i&i.ita) o o 
/(tfai.tfu) 

/ (2/31, J/32) 



< /(ai,£2) 




i = o 
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for decompositions yn + y 2 i + ?/3i = x\ and 2/12 + 2/22 + 2/32 = x 2 . This is so be- 
cause the solutions to the equation \t\ = t\ + t 2 = 0(mod3) are the multi-indices 
(1,2), (2, 1), (3, 3). Finally, if we set k = 3 and I = 2 the inequality (*) reads 

/ /(yn, 2/12, 2/23) \ 

f(yu, 2/22, 2/13) 

/(2/21, 2/12, 2/13) 

\ /(2/21, 2/22, 2/23) / 

< f(xi,x 2 ,x 3 )L 4: j = 

for decompositions 2/11 + 2/21 = £1, 2/12 + 2/22 = ^2 and j/13 + y 23 = £3- This is so 
because the solutions to the equation \t\ = t\ + t 2 + t% = (mod2) are the multi- 
indices (1, 1,2), (1, 2, 1), (2, 1, 1) and (2, 2, 2). 

Theorem 2.6 Let f :]0, ai[x • • • x]0, au[-^ R be a continuous, real function. If the 
function g : ]0, ai[x • • • x]0, afe[— > R defined by 

g(ri, ...,r k )=r^ 1 - ■ -r^ 1 f(n, ...,r k ) 

is matrix monotone of some index and order (I, . . . , I), then f is convex. 

Proof: We consider the simple root [3 = e 27Tl / 1 of the polynomial X 1 — 1 and set 

U = diag(/3P); =1 

which is a unitary matrix acting on C'. We introduce projections 

p = (u*yp u i 3 = 1,..., 1 

where P defined by 



p = \ 



is a one-dimensional projection acting on C l . We notice that 

1 



and consequently 



^ V =1 J P .o=l 



where E\ is the I x I identity matrix. The projections Pi, . . . , P/ are thus mutually 
orthogonal. Let xu, . . . , xu be real numbers in ]0, on{ and set 



Ci = diag^i^ ^ 
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The fc-tuple ((1 + le)x\, . . . , (1 + le)xk) is for a sufficiently small e > in the domain 
of g. We define 

Vji = x i /2 ( p i + e ) x \ /2 j = 1, . . .1; i = 1, . . . ,k 

and calculate 

yu H \-yu = x 1 / 2 (^2^ + x ^ 12 = ( 1 + le ) x * i = i,...,k. 

Since g is matrix monotone of index (Z, j) and order (I, . . . ,1) it follows that 
diag(g(y tl i,...,y tk k)) i , , J < #((1 + • • • , (1 + le)x k )L lk -i 

\ / |t|=J (modi) 

or inserting #(ri, . . . ,r k ) = rf 1 • ■■r^ 1 f(r 1 , ...,r k ) that 

diag^t/" 1 / 2 ® • • • O y^l 2 )f{y Sl i, ■ • ■ , y^/cXy" 1 / 2 <8> • • • <8> y" 1 / 2 )) , , . , jn 

\ / | s I = j (modi) 

< Ce k i( x i 1/2 ® ■ ■ ■ ® z~ 1/2 )/(c £ a;i, . . . , c £ x fc )(a^ 1/2 ® ■ ■ ■ ® x~ 1/2 )) 

\ / |t| = |s| = j (mod t) 

where c e = 1 + Ze. Multiplying to the left and to the right with the self-adjoint matrix 



diag^Y 2 <8> • • • <8> y]'X) . 



i (mod I) 

we obtain 



|s|— j (mod I) 



diag(/(y Sl i, . . . ,y Sfc fc)) 

< c ^ fe ((yl!i x i 1/2 ® • • • ® yl[ k x k 1/2 )f{ce x i, ■ • ■ , c £ x fe ) x 

We introduce for s» = 1, . . . , Z and i = 1, . . . , fc the Z x Z matrix 

Q Si i = (^f^V^V • 

ihI-IihI 1 " 92 /p,9=l 
It is an easy calculation to show that Q Sii is a projection and that 

X i — ' J ~Qsii s i — 1, . . . , Z, 2 1, . . . , h. 

Multiplying the above inequality from the left and the right with the projection 
diag(<3 Sl i <g> • • • <g> Q Sk k) 

V / \s\—j (modi) 
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and letting e tend to zero we thus obtain 
'x n -\ h% x lk H h% 



diag^Qsa 



Qs k k 



( xik H h ffzfc 

"A * 

((x} /2 P tl ® • •• ® xl /2 P tk )f(x u . . .,x k )(P sl x\ 
where we used that 



(4) < f x\i j ha^ii 

- V I 

,1/2 p. - 1/2 



1/2 



^ fc 4 /2 ) 



=j (mod I) 



\t\ — \s\—j (mod I) 



1/2 -1/2 
/M X i 



i'lt 



1/2 



QuiX r 



1/2 _ / glj j \-Xli 



I 



-1/2 



x \ /2 Pu 



as e tends to zero. We notice that (|J) is an l k 1 x Z fc 1 block matrix inequality of 
l k x Z fc matrices. Let us in order to examine the inequality calculate the entry 



4 /2 Pu 





}i /2 p t :_ 




x k r 4 




pq 




piqi 




Pkqk 



for p = (pi, . . . ,p k ) and q = (qi, ...,q k ) with pi,... ,Pfc,gi, . • . , qk = 1, We 
proceed to calculate the entry 



(a;J /2 F tl ® ••• ® xl /2 P tk )f(xii . . .,x k )(P Sl x[ 



1/2 



= J2 i- k p (u - pyt 4(1 ■ ■ ■ 4(1 ■ ■ ■ > x k){PsA 12 ® • • • ® PsA /2 ) 



u 1 ,...,u k = l 

I 



' P"' P ^ X p{l ' ■ ' X pikf ( X uil> ■ ■ ■ 1 x u k k)l •* X q[l''' X J h k 



ui,...,u k =l 



„.* 1/2 1/2 1/2 



,V2 



] — 2k l oq-s—p-i -,- - - - 

4 P X pil X Pk k X qil X q k k 



(3 {t - s) - u f(x ui i, 



i X Uk k) 



Ml,...,ltfe=l 



where we used that f(x%, . . . , is a diagonal matrix with f(x Ul \, . . . , x Ukk ) as the 
itth diagonal entry, and finally calculate the diagonal entry 



91 1 



We obtain from (Q) an inequality between / fe_1 x matrices by retaining the (t, s)- 
entry in each (t, s)-block on both sides of the inequality and discarding all other 
entries. We then insert the entries calculated above in the inequality obtained in this 
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way and get 



/( ; >•••, ; J [[(xn + ■ ■ ■ + x u ) L x 

^ ' i=l 

diag(x si i ■ • -x Skk 

k 

< r k Y[(x u + --- + x li )- 1 x 



| s I —j (mod I) 



2 = 1 



E /3 s - s - t -* +(t - s) -"/K 1 i ) ---,^ fe )) 

,...,« fe =i / 



/ 1/2 1/2 1/2 1/2 

I X t!l ' ' ' x t k k X s 1 l ' ' ' X s k k 

\ u u ...,u k = i / \ t \ = \s\=j (modi) 

Multiplying from the left and from the right with the self-adjoint matrix 
diag (x- 1 / 2 ■ ■ ■ x-H 2 f[(x u + ■■■+ x H yA 

V »=1 / |s|=j (modi) 

we obtain 

, ( x\\ H h in a;ife H \-xi k \ „ 

H z i J^- 1 

(5) 

«!,■■■, «fc=l 



|t| = |s|=j (mod I) 



We define for each u = (iti, . . . , Ufc) with . . . , Uk = 1, . . . , Z an l k 1 x l k 1 matrix 
II U by setting 

(6) n u = r^-v (ps-s-t-t+(t-s>A 

\ I \t\ = \s\=j (modi) 

It is an easy calculation to show that the matrices Tl u are self-adjoint projections, 
and the inequality (^|) can in terms of these projections be written as 

, f xu H \-xn xik H hx Zfc \ 

f { 1 1 J^- 1 

(7) j l 

— y ^ ^ f(XtLili • • ■ ) ^u fc /e)ll« ■ 



Because of 



it follows that 



j2 ^c«-*)-« = i*j t 



ui,. ..,«*=! 



(8) £ n u = 

Ul,...,«fe=l 
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Since each index (t, s) in each II n satisfy \t\ = \s\ = j (modi), confer equation (Jg), it 
follows that II„ = lit, for each v on the form 

(9) v = (vx, ... ,v k ) = (u± + i (modi), ... ,u k + i (modi)) i = 0, 1, . . . , I — 1. 

We also notice that for each u there are exactly I different indices in ([)]) . It follows that 
each projection is counted I times in the sum (^|). Two projections are consequently 
either orthogonal, or identical with their indices connected as in (^|). Setting u = 
(1, . . . , 1) and multiplying (0) with n u we obtain 

/ a:ii-| 1- xn xi k -\ 1- Xlk \ 

I I ) Uu 

< j(f(xxx, . . . ,xik) H h f(xn, . . . ,xik)jTlu- 

Therefore / is convex. QED 

A matrix monotone function may tend to minus infinity as the argument of the 
function approaches a point located on an axis, but it cannot go too fast. 

Corollary 2.7 Let g :]0, ai[x •••x]0,afe[— > R be a continuous function, which is 
matrix monotone of some index (l,j) and order (I, . . . , I). To each subset of the domain 
of g of the form ]0, /?i[x • • • x]0, f3k[ where (3\, . . . , (3k < oo, there is a constant C > 
such that 

C 

g(n,...,r k ) > (n,...,r k ) e]0,/?i[x ■■■x]0,/3 k [. 

Proof: The function / :]0, ai[x • • • x]0, a k [-^ R given by 

f(n, ■■■,r k )=r 1 - --rkg^x, ...,r k ) 

is convex by the preceeding theorem, and it is therefore bounded from below on 
bounded subsets of the domain. QED 

To proceed, we need the following slight generalization of |H, Theorem 1.2]. 

Theorem 2.8 Let f be a real, continuous function of k variables defined on the 
domain I\ X • • • X I k where I\, . . . ,I k are intervals containing zero and let (l,j) be 
any index. The following statements are equivalent: 

(i) f is operator convex and f(r\, . . . , r k ) < if ri = for some i — 1, . . . , k. 

(ii) The operator inequality 

diag(f(a* Sll xia Sll , a* Skk x k a Skk )j 
((at ll ®---®at k )f(x 1 ,...,x k )(a ail ®---®a Sh k)) 

\ * / \t\ = \s\=j (modi) 



\s\—j (mod I) 



< 



is valid for all unitary rows = (an, . . . , an) of length I acting on any Hilbert 
space Hi for i — 1, . . . , k and all k-tuples (x\, . . . ,x k ) of self-adjoint operators 
in the domain of f acting on H\ , . . . , H k . 
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(in) The operator inequality 

diag(f(p Sl ixip Sl u . . . ,p Skk x k p Skk ) 

\ / \s\—j (mod 1) 

< ( (Pta ® • • -®Pt k k)f{xi, ■ ■ . ,Xfc)(p ai i <S> • • • ®p Sfc fc)) 

\ / — |s|— j (mod /) 

is va/id /or aZ/ partitions of unity pn + ■ ■ ■ + Pu — 1 or* any Hilbert space Hi 
by orthogonal projections for each i = 1, . . . ,k and all k-tuples (xi, . . . , x k ) of 
self-adjoint operators in the domain of f acting on Hi, . . . , H k . 

The indices s, t in (ii) and (iii) are multi-indices of the form s = (si, . . . , s k ), where 
Si = 1, ... ,1 for i = 1, . . . , k with weight \s\ = si + • • • + sit- 
In the reference m the sufficiency of (it) and (iii) in order to obtain (i) were 
only established for indices of the form (1,0). However, rewriting of the original proof 
shows, mutatis mutandis, that the inequalities are indeed sufficient for the operator 
convexity of / for any index. The theorem above is stated for more general domains 
of the function / than in the original reference, cf. the discussion in the survey 
article || . It has the following version for functions of matrices S . 

Theorem 2.9 Let f be a real, continuous function of k variables defined on the 
domain I\ x ■ • • x I k where I\, . . . , I k are intervals containing zero and let (l,j) be any 
index. Let (m, . . . ,n k ) be a k-tuple of natural numbers and consider the statements: 

(i) f is matrix convex of order (Ini, . . . , ln k ) and f(r\, . . . , r k ) < if ri = for 
some i = 1, . . . , k. 

(ii) The matrix inequality 

diag(f(a* Sl iXia Sl i,...,a* k x k a Skk )) 

\ "J |s|=j (modi) 

< ((ali®---®at k )f(x u ...,x k )(a Sl i®---®a Skk )) 

\ / \t\ — \s\—j (modi) 



is valid for all unitary rows ai — (an, . . . , an) of length I acting on a Hilbert 
space Hi of dimension rii for i = 1, . . . , k and all k-tuples {x\, . . . , x k ) of self- 
adjoint operators in the domain of f acting on Hi, . . . , H k . 

(iii) The matrix inequality 

diag(f(p si ixip si i, . . . ,p Skk x k p Skk ) 

< [ (Ptil ® ' ' -®Pt k k)f(xi, ■ ■ .,X k )(p Sl l ® ■ ■ ■®Ps kk )) 

\ / |t| = |s| = j (modi) 

is valid for all partitions of unity pn + ■ ■ ■ + pu = 1 on a Hilbert space Hi of 
dimension Hi by orthogonal projections for each i = 1, . . . , k and all k-tuples 
(xi, . . . , x k ) of self-adjoint operators in the domain of f acting on Hi, . . . , H k . 



| s I —j (mod I) 
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(iv) The matrix inequality 

diag(f{p Sl ixips 1 i,...,Ps k kXkPs k k)) i , , 

\ / \s\—j {modi) 

< ((Pt 1 i®---®Pt k k)f(xi,...,x k )(p Sll ®---<S>p Skk )) 

\ / \t\ — \s\—j {modi) 

is valid for all partitions of unity pn + • • • + pn = 1 on a Hilbert space Hi of 
dimension In^ by orthogonal projections for each i = 1, . . . , k and all k-tuples 
(xi, . . . , Xf.) of self-adjoint operators in the domain of f acting on Hi, . . . , H k . 

(v) f is matrix convex of order (ni, . . . , n k ) and f(ri, . . . , r k ) < if Ti — for some 
i = 1, . . . , k. 

The implications (i) (ii) =>• (iii) and (iv) =>■ (v) are then valid. 

The indices s,t in (ii), (iii) and (iv) are multi-indices of the form s — (si, . . . , Sk), 
where Si = 1, . . . , I for i — 1, . . . , k with weight \s\ = Si + • • • + s&. Since the Hilbert 
spaces in (ii) are finite dimensional, it follows that any row = (an, ■ ■ ■ , an) satis- 
fying condition (|3|) is unitary. 

Theorem 2.10 Let f : [0, ai[x • • • x [0,afc[— > R be a continuous real function such 
that f(ri, . . . , rfc) < if Ti — /or some i = 1, . . . , k. The constants ol\, . . . , may 
be plus infinity. If f is matrix convex of order (Ini, . . . , Ink) for some integer I > 2 
and some k-tuple of natural numbers (ni, . . . , n k ), then the function 

g(ri, . . . ,r k ) = rf 1 • ■■r^ 1 f(ri,.. .,r k ) (n,.. . ,r k ) e]0,ai[x • • • x]0,a fc [ 

is matrix monotone of index (I, j) and order (ni, . . . , n k ) for j — 0, 1, ... ,1 — 1. 

Proof: Let (xi, . . . ,x k ) be any fc-tuple of positive invertible operators in the domain 
of / acting on Hilbert spaces Hi, . . . , H k of dimensions ni, . . . ,n k and let 

yu H 1- yu = Xi 

be any decomposition of Xi of length I for each i = 1, . . . , k. We set 

_ -1/2 1/2 _ 1 , . _ 1 , 

a sii — Xi y Sii Si — 1, . . . , f, I — I, . . . , K 

and observe that 

Vsii O Si iXiO Si i Si — \ , . . . ,1 , i — \ , . . . , k . 

If / is matrix convex of order (Ini, . . . , ln k ) we may apply Jensen's matrix inequality 
for functions of several variables, cf. Theorem |2.9| (») (ii). For each j = 0, 1, . . . , I— 1 
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we have 

diag(/(y si i, . . -,y Skk )) 
= diag(^/(a* ll a;ia Sl i, . . . , a* Skk x k a Skk )j 



\s\—j (mod I) 

)) 

\s\—j (mod I) 



< («i ® • • • ® a* tkk )f(xi, x k ){a Sll <g> ■ ■ ■ <g) a Skk )) 
= [ivlii^i 112 ® • • • ® yllk x k 1/2 )f(xi, ■■■,Xk)x 

= ((vlS ® • • • ® ylS)9(xi, ■ ■ • , ^X^ 2 ® • • • ® 2/^)) 



|£| = |s| = 7 (modi) 



\t\=\s\=j (modi) 



\s\—j (modi) 



\t\ — \s\—j (mod I) 

and multiplying to the left and to the right with the self-adjoint operator 

C = diag(y- 1 / 2 ®---®y- 1 / 2 ] 
in the above inequality, we obtain 

diag((2/ s " 1 1 1 (8>---®y7 fc 1 fc )/(y Sl i,...,y Sfc fe)) < (3(2:1, • • • ,»*)) . ., .,, 

V / |s|=j (modi) \ / \t\ = \s\=j (modi) 

or equivalently 

diag(g(y Sl i,...,y Sfcfc )) 1 , < 9(^1, ■ ■ ■ , x^Lp-i 

\ / I s I —j (modi) 

showing that g is matrix monotone of index (l,j) and order (m, . . . , rife). QED 

Theorem 2.11 Le£ / : [0, ai[x • • • x [0, a k [— > R be a continuous, real function and 
suppose the function 

0(ri, ...,7-fe) = r^ 1 • ■■r^ 1 f(ri, ...,r k ) (n, . . . ,r fc ) e]0,ai[x • • • x]0,a fc [ 

is matrix monotone of some index and order (lm, . . . , ln k ). Then the following 
statements are valid: 

(i) /(ri, . .. ,rfe) < ifri = /or some i = 1, .. ., k. 
(ii) f is matrix convex of order (m, . . . , rife). 



Proof: Since g is an increasing function in each coordinate, cf. Proposition 2.5, the 
first statement follows. 

Let {x\, ... ,x k ) be a fc-tuple of positive invertible operators in the domain of / 
acting on Hilbert spaces Hi, . . . , H k of dimensions Ini, . . . , ln k and let 

Pli H \-pu = l i = 1,.. . ,k 
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be resolutions of the identity on Hi of length I. We choose a positive e such that 
(1 + le)x is in the domain of / and set 

Usii — %J (Psii ~t~ &i = 1, . . . , l] % = 1, . . . , k. 

We consider the decompositions 

yu-\ h yu = (1 + le)xi i = l,...,k 

and use the assumption to obtain 

diag (g{y Sl i,... ) < g((l + le)xi, . . . , (1 + le)x k )L lk -i. 

V / I s I —j (modi) 

We introduce the diagonal block matrix 

C = d\a,g[x\ /2 {p Sl i + e)xip Sl i ® • • • ® x\ /2 {p Sk k + s)x k p Skk ) 

\ s \— 3 {modi) 

and multiply to the left with C* and to the right with C in the above inequality to 
obtain 

diag((p 81 i:ci(Paii + e)x\ /2 <g> • • • <g> p Skk x k (p Sk k + e)x]! 2 )g{y si i, ■ ■ -,y Sk k)x 
(x 1 1 /2 (p Sl i+e)x 1 p Sll ®---®xl /2 (p Skk +e)x kPskk )) 

/ \s\—j {mod I) 

< ((Ptii^ibtii +e)x\ /2 <g> • • • <g> pt k kX k (p tk k + e)xl /2 )g({l + le)xi, . . . , (1 + Ze)x fc ) 
x(^i /2 (^ii + e)ziP Sl i ® • •• ® 4 /2 (Ps fc fc + e)zfcP Sfe fc)), 

/ |t| = |s|=j (mod i) 

Inserting 

■ • ■ .l/sfcfe) = (^r 1/2 b S ii + er^i" 172 ® ■ ■ ■ ® x k 1/2 (Ps k k + e)- 1 x~ 1/2 )x 

f(x\ ,2 {p Sll + e)x\ /2 1 . . . ,xl /2 (p Skk + e)xl /2 ) 

and 

g((l + le)xi, . . . , (1 + le)x k ) = 

(1 + le)- k (x- 1/2 ® ■ ■ ■ ® x; 1/2 )/((l + ie)x!, . . . , (1 + le)x k )(x~ 1/2 ® • • • ® x~ 1/2 ) 
in the inequality, and then letting e tend to zero we obtain 

diag((p Sl ia;} /2 <g> • •• ®p Sfe fca;fe /2 )/(a;i /2 p Sl i:Ei /2 , . . . , x]/ 2 p Skk x]/ 2 )x 
(x 1 1 /2 p Sll x lPsil (g>---(g)x 1 k /2 p Skk x k p Skk )) 

/ |s|— j (modi) 

< {(Pt 1 ixipt 1 i ® • • ■®Pt k kX k p tkk )f(x 1 , . . .,x k )x 

[p Sl lXlp Sl l <g> • • • <g> Ps k kXkPs k k)) 

' \t\ = \ s \=3 (mod V) 
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| s I —j (mod I) 



\t\ — \s\ —j (modi) 



The identity 

p / 1/2 1/2 1/2 l/2w 1/2 _ „ 1/2 x 

/(V Psii^i Psfcfc^fc Hx/ p ai i®---®»fe Psfefcj = 

Ol /2 Psil <S> • • • ® xl /2 p Sh k)f(p Sl lXlp Sl l, . . . ,Ps k kX k p Sk k) 

follows by first considering polynomials and then applying Weierstrass' approximation 
theorem. Inserting the identity in the inequality above we obtain 

diag((p Sl iiz;ip Sl i ® • • ■ <g> p Sk kXkPs k k)f{p Sl iXip sl i, ■ ■ ■ ,Ps k kX k p Sk k)x 

(p sl iXiP sl i <8> • • • ®Ps k kXkP Sk k)^j 
< {{ptiixiPtti ® ■ ■ ■ <8>Pt k kXkPt h k)f(xi, ■ ■ -,Xk)x 

(p Sl lXip Sl i (g) • • • <S> Ps k kXkPs k k)) 

and hence 

diag((p si i ® • • ■®p sk k)f(p si ixip si i, ■ ■ ■ ,p Sk kXkPs k k)(p sl i ® • • • ®P Sk k)j 
< ((Pt 1 i®---®Pt k k)f(xi,...,x k )(p Sll <g)---<g>p Sk k)) 

\ / |i|^|s|— j (modi) 

Because of (i) we obtain 

f{Ps 1 lXlPs 1 l, . . ■ ,p Skk XkPs k k) 

< (P S1 1 <8> • • • ®Ps k k)f{p Sl \XlP Sl l 1 . . . ,Ps k kXkPs k k){p Sl l <8) • • • ®Ps k k) 

and consequently 

di&g(f(Ps 1 ix 1 p Sll ,...,p Sk kXkPs k k)) ti r jn 

\ / \s\—2 [mod I) 

< {{Phi <8 • • -®Pt k k)f(xi,. . .,x k )(p Sl i ® ••■ <8iPs k k) ) 

\ / \t\=\s\=j (modi) 



\s\—j (mod I) 



which is Jensen's matrix inequality. We thus deduce, cf. Theorem 2.E (iv) => (w), 
that / is matrix convex of order (m, . . . , n^.). QED 

One may think that the proceeding theorem, which ensures matrix convexity of 



/, could replace Theorem 2.6 which with similar conditions only imparts ordinary 



convexity on /. However, it is essential in the proof of the preceeding theorem that / 



is defined also on the axes, while this is not required in Theorem 2.6. This problem can 
easily be overcome for functions of only one variable by making a small translation of 
the matrix monotone function g. This remedy is not available for functions of several 
variables, since the translation of the decomposition of an operator no longer is a 
decomposition of the translated operator, cf. equation (||). 

Corollary 2.12 Let g :]0, a\[x ■ ■ ■ x]0, a/c[— > R be a continuous real function. If g 
is matrix monotone of some index (l,j) and order (Irani, . . . ,lmn k ) for a natural 
number m and a k-tuple of natural numbers (n\, . . . , n k ), then it is matrix monotone 
of index {in, h) and order (m, . . . , nk) for h = 0, 1, . . . , m — 1. 
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Proof: The real and continuous function / defined by 

f(rx,...,r k ) =rx---r k g{rx,...,rk) < < a* for i = l,...,k 

is convex by Theorem [2.6| . It therefore extends to a continuous function 

/ : [0,ai[x • • • x [0,a fc [-> R, 

and it follows that f(ri, . . . , r k ) < if = for some i = 1, . . . , k. The function / 

i g is thus 
,m—l by 



is matrix convex of order (mni, . . . ,mrik) by Theorem 2.11. The function g is thus 
matrix monotone of index (m,h) and order (n\, 
Theorem WM. QED 



, rzfe) for h = 0, 1, 



Corollary 2.13 Let g :]0, oji[x • • • x]0, * R &e a continuous function. If g is 
operator monotone of some index, then it is operator monotone of all indices. 



Definition 2.14 We say that a continuous function g :]0, a%[x ■ • • x]0, ctk[-^ R is 
operator monotone, if it is operator monotone of some and hence operator monotone 
of all indices. 



Pro of (o f Theorem 1.2) : Th e statement follows by combining Theorem 2.10 , Theo- 
rem |2.1l| and Definition [2.14|. QED 



The simplest ex amp le of operator convex functions satisfying the boundary con- 
ditions in Theorem |l.2j are the negative constants. The function 

g(n,...,r k ) = -r^ 1 ■ ■ ■ rZ 1 



defined in the first (open) quadrant is thus operator monotone, cf . also Corollary [2/? . 
The set of operator monotone functions defined on a given domain is a weakly closed 
convex cone, but the constant function g(r\, . . . ,rfc) = 1 is not operator monotone 
for k > 2. This must indeed be so since the function (n, r2) — > r\r2 is not convex. 



Frank Hansen: Institute of Economics, University of Copenhagen, Studiestraede 6, DK-1455 
Copenhagen K, Denmark. 
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